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An approximation that is often used in fits to reactor and atmospheric 
neutrino data and in some studies of future neutrino oscihation experiments 
is to assume one dominant scale, Am^, of neutrino mass squared differences, 
in particular, Am^^^ ~ 3 X 10-3 eV^. Here we investigate the corrections 
to this approximation arising from the quantity Am'j.^i relevant for solar neu- 
trino oscillations, assuming the large mixing angle solution. We show that for 
values of sin (26I13) ~ 10^2 (ij^ ^j^g 

range of interest for long-baseline neutrino 
oscillation experiments with either intense conventional neutrino beams such 
as JHF-SuperK or a possible future neutrino factory) and for Am^^; ~ 10^^ 
eV'^, the contributions to — > oscillations from both CP-conserving and 
CP-violating terms involving sin^(Am^^;L/(4£^)) can be comparable to the 
terms involving sin^(Am^j^L/(4£')) retained in the one-Am^ approximation. 
Accordingly, we emphasize the importance of performing a full three-flavor, 
two-Am^ analysis of the data on fe oscillations in a conventional-beam 

experiment and ^ f^, i^e ^ oscillations at a neutrino factory. We also 
discuss a generalized analysis method for the KamLAND reactor experiment, 
and note how the information from this experiment can be used to facilitate 
the analysis of the subsequent data on i/^ oscillations. Finally, we con- 

sider the analysis of atmospheric neutrino data and present calculations of 
matter effects in a three-flavor, two-Am-^ framework relevant to this data and 
to neutrino factory measurements. 
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I. INTRODUCTION 



There is increasingly strong evidence for neutrino oscillations, and thus neutrino masses 
and lepton mixing. All solar neutrino experiments that have reported results (Homestake, 
Kamiokande, SuperKamiokande, SAGE and GALLEX/GNO, SNO) show a significant deficit 
in the neutrino fluxes coming from the Sun [1]. This deflcit can be explained by oscillations of 
the into other weak eigenstate(s). The currently favored region of parameters to flt this 
data is the solution characterized by a neutrino mass squared difi^erence 2 x 10~^ ^ ^^Lz ^ 
2 X 10^"^ eV^ and a large mixing angle (LMA), tan^ 6'i2 ~ 0.4 [1]- [3]. Solutions yielding 
lower-likelihood fits to the data include the small-mixing angle (SMA) solution, with strong 
Mikheev-Smirnov-Wolfenstein (MSW) matter enhancement [4], and ^rn^soi ~ O-^ ^ 
eV^, tan^ ^12 ~ 4 x 10~^, the LOW solution with Am^^; ~ 10~^ eV^ and essentially maximal 
mixing. 

Another piece of evidence for neutrino oscillations is the atmospheric neutrino anomaly, 
observed by Kamiokande [5], 1MB [6], Soudan [7], SuperKamiokande (also denoted SuperK, 
SK) with the highest statistics [8], and MACRO [9]. The SuperK experiment has fit its 
data by the hypothesis of — > Vt oscillations with Am^^^ ~ 3 x 10~^ eV^ and maximal 
mixing, siv?26atm = 1- The possibility of Ug oscillations involving light electroweak- 

singlet ("sterile") neutrinos has been disfavored by SuperK, and the possibility that — > Ve 
oscillations might play a dominant role in the atmospheric neutrino data has been excluded 
both by SuperK and, for the above value of Am^^^, by the Chooz and Palo Verde reactor 
antineutrino experiments [10,11]. The K2K long-baseline neutrino experiment between KEK 
and Kamioka has also reported results [12] which are consistent with the SuperK fit to its 
atmospheric neutrino data. The LSND experiment has reported evidence for z/^ — > Vf, and 

~^ oscillations with Am^^^^) ~ 0.1 — 1 eV^ and a range of possible mixing angles [13,15]. 
This result is not confirmed, but also not completely ruled out, by a similar experiment, 
KARMEN [14]. The solar and atmospheric data can be fit in the context of three-fiavor 
neutrino mixing, and we shall work within this context. 

The fact that these inferred values of neutrino mass squared differences satisfy the hier- 
archy lAm^Q^I << |Am^j^| has led to a commonly used approximation in fits to the reactor 
and atmospheric neutrino data, and in studies of CP-conserving effects in terrestrial neu- 
trino oscillation experiments. In this approximation, which we shall denote the one-Am^ 
approximation (IDA), one neglects Am^o^ compared with Am^^^. For certain neutrino os- 
cillation transitions, such as Ur, this is an excellent approximation. It is worthwhile, 
however, to have a quantitative cvahiation of the corrections to this approximation and a 
determination of the ranges of parameters where these corrections could become significant. 
Indeed, for sufficiently small values of the lepton mixing angle 9is (defined below in eq. 
(2.2)), e.g., sin^2^^i3 ~ 10"^, and sufficiently large values of ^rn^^i, e.g., Am^^^ ~ 10~^ eV^, 
this approximation is not reliable for certain oscillation channels such as z/g. Here we 
are referring to CP-conserving quantities; the one-Am^ approximation is, of course, not used 
for calculating CP- violating quantities since neglecting Am^^^ is equivalent to setting two of 
the neutrino masses equal (in the standard three-active-fiavor context), which allows one to 
rotate away the CP-violating phase that would appear in neutrino oscillation experiments 
and hence renders CP- violating quantities trivially zero. Since the values of sin^(2^i3) and 
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Am^^; for which the one-Am^ approximation breaks down are in the range of interest for 
future experimental searches for i/^ — > i/^ via both conventional neutrino beams generated by 
pion decay and via neutrino beams from neutrino "factories" based on muon storage rings, 
this complicates the analysis of the sensitivity and data analysis from these experiments. 
Many fits have been performed of solar and atmospheric data, e.g. [2,3], [16]- [19]. One 
salient result is that in Ref. [18], the usual fit to the SuperK atmospheric neutrino data with 
a single Am^^^ ~ 3 x 10~^ eV^ and siii^ {26atm) = 1 is compared with a very different fit 
with two equal mass squared differences, Ar7i32 = Am|j = 0.7 x 10~^ eV^, and it is argued 
that although the for the latter is worse, it is still an acceptable fit. This suggests that 
one should carefully assess corrections to the one-Am^ approximation in studies of neutrino 
oscillations. 



II. GENERALITIES ON NEUTRINO MIXING AND OSCILLATIONS 

A. Mixing Matrix and Oscillation Probabilities 

In this section we briefiy record some standard formulas for neutrino oscillations that we 
shall use. In the framework of three active neutrinos, the unitary transformation relating 
the mass eigenstates Ui, i — 1,2,3, to the weak eigenstates Ua is given by Ua — J2^=i Uai^i 
where the lepton mixing matrix is 

U = R2^KR^^K*RnK' = (2.1) 
^ C12C13 S12C13 Si3e '^^ \ 

= -S12C23 - ci2S23Si3e*'' C12C23 - si2S23Si3e"^ S23C13 \k' (2.2) 

V S12S23 - Ci2C23Si3e"^ -C12S23 - Sl2C23Sl3e"^ C23C13 / 

Here Rij is the rotation matrix in the ij subspace, Cij — cos%, — sin%, K — 
diag{e'^^, 1, 1) and K' — diag{l, e"'^, e"'^) involves further possible phases (due to Majorana 
mass terms) that do not contribute to neutrino oscillations (as can be seen from the invari- 
ance of the quantity Kah,ij below under neutrino field rephasings). One can take O-ij G [0, 7r/2] 
with 5 e [0, 27r). The rephasing- invariant measure of CP violation relevant to neutrino oscil- 
lations is given by the Jarlskog invariant [20] determined via the product Im{UijUknU*nU^j), 



- sin(2^i2) sin(2^23) sin(2^i3) cos ^13 sin 5 (2.3) 

8 

In vacuum, the probability that a weak neutrino eigenstatc Va becomes Vh after propa- 
gating a distance L (assuming that E » m{ui) and the propagation of the mass eigenstates 
is coherent) is 

3 

Pil^a T^b) = 5a& - 4 ^ Re{Kab,ij) sin^ 

i>j=l 

3 

+ 4:J2 Im{Kab,ij) sin (f)ij COS (f)ij (2.4) 
i>j=i 
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where 

Ka,,ij = Ua.U;,U:^U,j (2.5) 

AmJ = m{uif - m{vjf (2.6) 

and 

*i = (2.7) 

We recall that for the CP-transformed reaction and the time-reversed reaction 

^ J^a, the oscillation probabilities are given by eq. (2.4) with the sign of the Im{Kab,ij) 
term reversed. Further, by CPT, P(i?ft — > Vg) — P^i'a ~^ ^b) so that, in particular, P(Pa ~^ 
Vg) = P{i^a ^a)- It is straightforward to substitute the elements of the lepton mixing 
matrix (2.2) and evaluate the general formula (2.4) for each of the relevant transitions. 
For the special case P{i'a ^a), eq. (2.4) simplifies to 

3 

Piiya^jya)^i-4: J2 \Ua^\^\UaJ\^Sm^4>iJ (2.8) 
i>j=l 

We recall the elementary identity 

Am32 + Amli + Aml^ = (2.9) 

so that in general, three- flavor vacuum oscillations depend on the four angles 612, O23, ^13, S 
and two Am^'s, which can be taken to be Am32 and Am2i. The currently favored regions 
for 621 and Avn^^ are determined primarily by the solar neutrino data. Here one can take 
Am|i > with 6*12 G [0,7r/2] [21]. To distinguish between the first and second octants, 
the parameter regions allowed by these fits to the solar data can be expressed in terms of 
Arr?2\ > and tan^6'2i. 

The commonly used one-Am^ approximation is then based on the hierarchy 

Amli = Ami, < l^^lil ~ l^m^^l = \AmlJ (2.10) 

However, as mentioned above, the solar data itself or in combination with atmospheric and 
reactor data allows for rather large values of Am'^^i, up to ~ 2 x 10""^ eV^ or even somewhat 
higher. 

B. Two-flavor Oscillations 

In the case of oscillations of two fiavors, the oscillation probability in vacuum is, in an 
obvious notation, 

P{i^a T^b) = sin^(2e) sin^ 0^ (2.11) 
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C. Three-Flavor Oscillations with One-Am^ Dominance 



With the hierarchy (2.10), one has the following approximate formulas for vacuum oscilla- 
tion probabilities relevant for experiments with reactor antineutrinos, atmospheric neutrinos 
and CP-conserving effects in terrestial long-baseline oscillation studies: 

P{v^ Ur) = 4|C/33HC/23|^sin2(/)„t„ 

= sin2(2^23) cos^ ^13 sin^ cf)^,^ (2.12) 

Pil'u. ^e) = ^\Ul2,\^\U23,\^ (t)atm 

= sin2(2^i3) sin' ^23 sin' 0atm (2.13) 

P{Ve I^t) = 4|f/33p|f/i3psin'0at„ 

= Sin2(2^13) cos' ^23 sin' (j)atm (2.14) 

Since this one-Am' approximation removes CP-violating terms, it also implies that the 
oscillation probabihties for the CP-transformed and T-reversed transitions are equal to the 
probability for the original transition, P{i?a h)iDA = Pin i^ajiDA = -P('^a n)iDA- 

For the analysis of data on reactor antineutrinos, i.e. tests of Pg ~^ ^e, the two-flavor 
mixing expression is P(t'e — >■ i>e) = 1 — P(t'e 7^ ^e); where 

Pil^e ^ ye) = Sin'(2^,eactor) siu' ( ^'"^'g"'^^ ) (2.15) 

where Aml^^^^^j, is the squared mass difference relevant for Pg D^. Combining (2.13) and 
(2.14), we have, in this approximation, 

Oreactor = ^13 , ^■ml^actor = ^"^32 (2-16) 

For the analysis of atmospheric data with the transition favored by the current data, 
letting 

P{v, vr) = sin'(2e„,^) sin'(^^^) (2.17) 

one has, using (2.12), 

sin'(2^„,„) = sin'(2^23) cos^ ^13 (2.18) 

and Am'j^ as given in (2.10). Since the best fit value in the SuperK experiment is 
siv?{26atm) = 1, it follows that 

^13 « 1 (2.19) 



4 



and hence 



Oatm - 023 (2.20) 

For the K2K experiment, using ^13 << 1, one has 

P{u^ ^ i/^) ~ 1 - sin2(2^23) sin^ (2.21) 

All of these vacuum oscillation probabilities are independent of the sign of Am^^^, just 
as in the two-flavor vacuum case. However, the symmetry 9 — > 7r/2 — 6 of the two- flavor 

vacuum case is no longer present. For ^13, one can immediately infer that this angle is near 
rather than near to tt/2 from the fit to the atmospheric data, as noted above. For 623, the 
transformation 623 — >■ 7r/2 — 623 leaves the expression for P(t'^ Vt) in (2.12) invariant and 
interchanges the values of the oscillation probabilities P{ve ~^ ^fj.) and P{i'e ~^ ^t)- Because 
we know that the value of 623 is close to 7r/4, this interchange does not make a large change 
in these probabilities P{i'e ^/i) and P(z/e — * i^r)- The atmospheric data places an upper 
bound on the transition P(z/^ i^e), and the fact that this is small is implied by the fact 
that 6*13 << 1, so that this upper bound does not determine how large the sin^6'23 factor in 
(2.13) is and hence whether $23 is slightly below or slightly above 7r/4. 



III. GENERALIZED ANALYSIS OF REACTOR ANTINEUTRINO DATA 

The general three-flavor, two-Am^ (i.e., two independent Am^) formula for antineutrino 
survival probability that is measured in reactor experiments such as Gosgen, Bugey, Chooz, 
Palo Verde, and KamLAND is P(z/e z/g) = 1 — P(t'e 7^ i^e); where, using (2.8), we have 
(note that by CPT, P(z/e 7^ u^) = P{u^ z/^)) 



.ej\ Sm (pij 



3 

i>3=l 

= sin^(26'i3) sin^ 6*12 sin^ ^32 sin^(26'i3) cos^ 6*12 sin^ ^31 

+ sin2(2^i2) cos^ ^13 sin^ 02i (3.1) 

Matter effects are negligible for these experiments. Let us consider the results from the Chooz 
experiment, since these place the most stringent constraints on the relevant parameters. This 
experiment obtained the result [10] 

R = ^^"'''^^'^ = 1.01 ± 0.028(stot.) ± Q.m{sys.) (3.2) 

calculated 

From this and the agreement between the measured and expected positron energy spectra, 
this experment set the limit 

P{v, -h Ve) < ech (3.3) 
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where ech — 0.05. Within the usual context of two-flavor mixing described by (2.15), the 
Chooz experiment reported the 90 % confidence level (CL) hmits 

sin^ 29reactor <0-^ for large Aml^^^^ (3.4) 

and 

|A^4actorl < 0.7 X 10-3 eV' for sin'' 29reactor = 1 (3.5) 

In the one-Am^ approximation, the first two terms of eq. (3.1) combine to make the term 
sin^ 2^13 sin^ 4>atm so that 

Sin^ 26'i3 = Sin^ 29reactor < 0.1 (3.6) 

Since each of the three terms Tj in the general equation (3.1) is positive-definite, we 

have, in an obvious notation, Tj < ech- Applying this to the third term and using the fact 
that ^^13 << 1, we obtain an upper bound on sin^ 29i2 sin^ 02i- From the plot of the Chooz 
excluded region [10], we infer the pair of bounds 

f 0.7 X 10-3 eV^ for sin^ 29i2 = 1 
Ami, < (3.7) 
[ 0.9 X 10-3 eV^ for sin^ 29u = 0.8 

where the second upper bound applies for the central value of sin^ 29i2 — 0.8 in the LMA 
solution, 

LMA(central) : tan^ 9i2 = 0.4 , Amj, = 0.5 x 10-^ eV^ (3.8) 

For the central values of the LMA solution in (3.8), using L — 1 km and a typical energy 
E A MeV, the third term in (3.1) has a value of about 1.3 x 10"^, which is neghgibly 
small. This increases to about 2 x 10-^ for Amli = 2 x 10"^ cV^, which is again negligible 
compared with the range of z/g disappearance, ~ 0.05, probed by Chooz. 

Next, we consider the KamLAND long-baseline reactor experiment, which will use a 
liquid scintillator detector in the Kamioka mine to measure P^p e+n events initiated by 
Pe's from a number of power reactors and thereby test the LMA solution to the solar neutrino 
deficit and is expected to begin data-taking in 2001 [22]. The power reactors are located 
at various distances from 140 km to 200 km from Kamioka. It has been estimated that, in 
the absence of oscillations, a total of 1075 UeP events per kton-yr will be recorded, 

and of these, 348 events per kton-yr will arise from the single most powerful reactor, the 
Kashiwazaki 24.6 GW (thermal) facility a distance L — 160 km away [22]. For the conditions 
of this experiment, \(f)3j\ » 1 for j = 1,2, so that the sin^ 03^ factors average to 1/2 over 
the Pg energy spectra from the reactors, and hence (3.1) reduces to 

P{l^e 7^ i^e)Kam. = \ Sm'(2^i3) + sin2(2^i2) COS^ ^13 siu^ 021 (3.9) 

The first term has a maximum value of 0.05. For the central LMA values in (3.8), the second 
term has a value of approximately 0.1 (almost independently of ^13, given the bound (3.6)). 
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Thus, if sin^(2^i2) and Am^^ are characterized by the LMA solution and if sin^(2^i3) is near 
to its current upper bound, then the two terms in eq. (3.9) would make contributions that 
differ only by about a factor of 2. One can thus distinguish several possible outcomes for 
the KamLAND experiment: 

• If this experiment sees a signal for oscillations of reactor PgS with P{Pe 7^ ^e) > 0.05, 
this implies that there is at least some contribution to this signal from the second term. 

• In general, if P{i^e 7^ ^e) is nonzero, then, from the overall deficiency in the rate, one 
would not be able to determine the relative contributions of each term in (3.9). Instead, 
for one pathlength, L, one would have to perform a three-parameter fit involving the 
parameters cj, j = 1,2,3, where Ci is a constant, representing the first term in (3.9), 
and Cj, j = 2,3 enter as C2 sin^(c3£'/(4L)), representing the second term in (3.9). Since 
the KamLAND detector is sensitive to Pg's from a number of different reactors at 
different distances, the actual fit to the data would be more complicated than this, 
but the oscillations would not, in general, be adequately described by a two-fiavor, 
one-Am^ formula, and cq. (3.9) would apply for the three-fiavor, two-Am^ analysis, 
given the size of Am^j^ inferred from the atmospheric data. Clearly, the actual results 
extracted from the data will depend on simulations that involve the determination and 
subtraction of reactor backgrounds [23] . 

• Finally, if the KamLAND experiment sees no signal for oscillations and sets the limit 
P(Pe Ue) < ^KL, then since each of the terms in (3.9) is positive-definite, one will have 
the bounds sin^(2^i3) < 2eKL and the usual excluded-region plot for the contribution of 
the second term. This will depend on the statistical uncertainties and the backgrounds 
and resultant systematic uncertainties, but, roughly speaking, it will have asymptotes 
sin^(26'i2) < 2eKL if 02i is assumed to be large enough so that sin^ 02i averages to 1/2 
over the reactor Pg energy spectra, and a corresponding bound of sin^ 02i < ^kl for 
maximal mixing, sin^ 2^12 = 1 (given that one knows that the cos^ ^13 factor is very 
close to unity) . 

IV. GENERALIZED ANALYSIS OF A LONG-BASELINE EXPERIMENT TO 

MEASURE ^ i/e 

In this section we shall discuss the general three-flavor, two- Am^ analysis of long-basehne 
accelerator experiments to measure — > using conventional beams and v^, — > or its 
conjugate using beams from a possible future neutrino factory based on a muon storage 
ring. There are several long-baseline accelerator experiments under construction to continue 
the study of neutrino oscillations after the pioneering work of the K2K experiment. These 
include the MINOS experiment from Fermilab to the Soudan mine, with L = 730 km, using 
a far detector of steel and scintillator and a neutrino flux peaked at £^ ~ 3 GeV [24]. In 
Europe, a program is underway to use a neutrino beam with ~ 20 GcV from CERN a 
distance L = 730 km to the Gran Sasso deep underground laboratory, involving the OPERA 
experiment and also plans for a liquid argon detector [25]. Third, the JHF-SuperK neutrino 
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oscillation experiment will use a i/^ beam from the 0.75 MW Japan Hadron Facility (JHF) 
High Intensity Proton Accelerator (HIPA) in Tokai, travelling a distance L — 295 km to 
Kamioka [26] . In a first stage, this would use SuperK as the far detector; a possibility that 
is discussed for a second stage involves an upgrade of JHF to 4 MW and the construction 
of a 1 Mton water Cherenkov far detector (denoted HypcrKamiokandc). The JHF-SuperK 
collaboration has stated that one of the goals of its first phase is to search for i/^ — >■ 
oscillations down to the level -P(^'^ — > Ve) ~ 0.003 by taking advantage of the excellent 
particle identification ability and energy resolution of SuperK for electrons and muons [26]. 
This is the sensitivity for a narrow-band beam with E = 0.1 GeV, which, assuming that 
|Am^2l = 3 X 10"^ eV^ maximizes the factor sin^032; the estimated sensitivity for a wide- 
band beam with E peaked at about 1.1 GeV is -P(t'/i I'e) ~ 10^^ [26]. This type of search 
will be pursued to some level also by the other long-baseline experiments. A number of other 
possible long-baseline z/^ Vf. oscillation experiments using intense conventional neutrino 
beams have been considered, with a variety of pathlengths [27]. 

A different approach that has been considered in detail is that of a neutrino "factory", 
in which one would obtain a very intense beam of and z/g's from the decays of /i~'s in 
a muon storage ring in the form of a racetrack or bowtie, and similarly a beam of and 
i^e's from stored /i"'"'s. These beams would have a definite and precisely understood fiavor 
composition and would make possible neutrino oscillation searches using very long-baselines 
of order 3000 km [28]- [33]. Typical design parameters arc E = 20 GcV for the stored /i"^ 
energy and L = 10^° ii decays per Snowmass year (10^ sec). With a stored /x~ beam, say, one 
would carry out a measurement of the survival probability via the charged current 

reaction yielding a final state and an appearance experiment with Pg yielding a 

final state yU"*", a so-called wrong-sign muon signature. It has been estimated that with a 
moderate-level neutrino factory, one could search for ^ti or its conjugate reaction down 

to the level sin^ (26*13) ~ 3x 10^^ [30]. For such long pathlengths, matter effects are important 
[33] and can be used to get information on the sign of Aml2- It may also be possible to 
measure leptonic CP violation using either a conventional beam or a beam from a neutrino 
factory. 

At the levels of sin^ 2^^i3 that will be probed, the one-Am^ approximation used in many 
planning studies may well be inadequate, and one should use a more general theoretical 
framework. The full expression for the ly^ — > oscillation probability in vacuum (matter 
effects are discussed below) is obtained in a straighforward manner from the formulas (2.4) 
and (2.2) and has the form, in a compact notation, 

P{j^n T^e) = 2sin(26li3)s23Ci3Si2(si2S23Si3 " CiaCaaC^) sin^ ^32 + 



2 Sin(26li3)s23Cl3Cl2(ci2S23Sl3 + S12C2ZC5) SVC? 031 - 
- 2 Sin(26'i2)c^3 Sl2Ci2(s?3S23 - C23) + Sl3S23C23(Sl2 " Ci2)c5 siu^ 021 

]- sin(26'i2) sin(26'i3) sin(26'23)ci3S5 sin ^32 cos ^32 - 



— sin cos 031 + sin 02i cos 02i 



(4.1) 



(As in (2.4) 



P(Py 



Ve) and P{i'e 



Vfj) are given by (4.1) with the sign of the sin (5 
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term reversed.) Now for sufficiently small Am|^, as would be true in the solar neutrino fits 
with SMA, LOW, or Am|i ~ 10~^ eV^ vacuum oscillations, and sufficiently large sin^(2^i3), 
subject to the constraint (3.6), the full eq. (4.1) reduces to (2.13), in which the oscillation is 
driven by the terms involving sin^ (patm = sin^ 032 — sin^ 03i. However, if Am2i and sin^ 2021 
are at the upper end of the LMA region, then the one-Am^ approximation can break down. 
As a numerical example, one can consider the parameter set sin^ 29i2 = 0.8, Am|j^ = 2 x 10^^ 
eV^, sin^ 2^13 = 0.01, 5 = 7r/6, with the usual central SuperK values Aml2 = 3 x 10~^ eV^ 
and sin^2^23 = 1- Further, take the JHF-SuperK pathlength L — 295 km and narrow-band- 
beam energy E = 0.7 GeV, and label this total set of parameters as set (a). Then, if one 
were to evaluate the z/^ oscillation probability using the one-Am^ approximation, again 
denoted IDA, eq. (2.13), one would obtain 

P{l^^, J^e) = 5.0 X 10"^ for set (a) with IDA (4.2) 

However, correctly including the contribution from the term involving sin^ 02i, using the full 
expression (4.1), one gets an oscillation probability that is more than twice as large as the 
one predicted by the one-Am^ approximation: 

P{u^ z/e) = 1.4 X xlO-2 for set (a) (4.3) 

This clearly shows that for experimentally allowed input parameters involving the LMA so- 
lar fit, and in particular, for a value of sin^ 2^13 that can be probed by the JHF-SuperK 
experiment and others that could achieve comparable sensitivity, the one-Am^ approxima- 
tion may not be valid. Thus, it is important that the KamLAND experiment will test the 
LMA and anticipates that, after about three years of running, it will be sensitive to the level 
Am^p; ^ 10~^ eV^ [22]. This information should therefore be available by the commissioning 
of JHF in 2007. The adequacy of the three-flavor theoretical framework will also be tested 
by the miniBOONE experiment within this period. If, indeed, the LMA parameter set is 
confirmed by KamLAND, then it may well be necessary to take into account three-flavor 
oscillations involving two independent Am^ values in the data analysis for the JHF-SuperK 
experiment and other z/^ — ^ 74 neutrino oscillation experiments that will achieve similar sen- 
sitivity. This point is thus certainly also true for long-baseline experiments with a neutrino 
factory measuring i/g — > i^^, Pe — P/* oscillations, since they anticipate sensitivity to values 
of sin^ 2^13 that are substantially smaller than the level to which the JHF-SuperK collabora- 
tion will be sensitive, and as one decreases ^^13 with other parameters held fixed, the sin^ 02i 
corrections to the one-Am^ approximation become relatively more important. 
In passing, we observe that in the limit 613 0, eq. (4.1) reduces to 

P{u^ Ue) = sin2(2^i2) cos^ 023 sm^ 021 for ^13 = (4.4) 

In this limit, the term involving sin^ 02i, rather than the terms involving sin^ 032 or sin^ 03i, 
are driving the i/^ — > Ug oscillations. 

V. 1^1^ ^ z/^ DISAPPEARANCE EXPERIMENTS 

All long-baseline accelerator neutrino experiments, including K2K, MINOS, CNGS, JHF- 
SuperK, and other possible ones such as CERN-Frejus and those that might involve UNO 
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and/or a neutrino factory, will perform a measurement of the survival probability. 

The one-Am^ approximation yields the result 

7^ ^^m) = 4(1 - it/^3r)it/M3rsin2 032 

= sin2(2^23)cos'^i3 + sin2(2^i3)sin4^23 sin^ (5.1) 

Since SuperK infers a maximal oscillation to fit its atmospheric neutrino data, and 

since this implies that ^^13 << 1, the second term in (5.1) is quite small compared to the 
first. As a numerical example, for sin^ 26'i3 = 0.01 and 6*23 = 7r/4, the ratio of the second 
to the first term in (5.1) is 2.5 x 10~^. The one-Am^ approximation is a very good one for 
this transition; for experiments such as MINOS and JHF-SuperK, the relative corrections 
are typically of order < O(10~^). 

VI. — 

This transition is more difficult to measure than z/^ since (a) the optimal neutrino 

energy to maximize the oscillation factor is below r threshold, and (b) even if this were not 
the case, the r is not observed directly. For completeness, however, it should be noted that 
again the term retained in the usual one-Am^ approximation, (2.14) may not be larger than 
the term that would describe this transition if ^13 = 0, namely 

P{y^ v^) = sin2(2^i2) sin^ ^23 sin^ 02i for ^13 = (6.1) 

This is the same as the expression for P(ve — > v^^ = P(z/^ — > z/g), eq. (4.4) under the same 
assumption, 6^13 = with the interchange of cos^6'23 and sin^6'23. 

VII. MATTER EFFECTS FOR NEUTRINO OSCILLATIONS WITH TWO 

RELEVANT Am^ SCALES 

In many experiments matter effects can be relevant. This is the case with solar neutrinos, 
atmospheric neutrinos, and future possibilities for 0(10"^) km baseline neutrino oscillation 
experiments using neutrino factories [30]- [38]. In these cases, oscillation probabilities are 
modified by the interaction of the neutrinos in the matter: i/^ and have the same forward 
scattering amplitude, via Z exchange, while Vg has a different forward scattering amplitude 
off of electrons, involving both Z and W exchange. This leads to a matter-induced oscillation 
effect when electron neutrinos are involved in the oscillations. 

In this case one needs to solve the evolution equation which includes the effects of the 
interactions with matter, which reads (for a generic two-generation case) 

^(''') = (^UM^U^ + V 

where 
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M 



Ve 





COS 9 sin 9 
— sin 9 cos 9 



(7.2) 



Here V" = Ve = V^GpNe where iVg is the electron number density and we have v^G^A^e 
[eV]= 7.6 X 10~^^FeP [g/cm^], where p is the mass density and Yg is the average electron 
fraction of the matter. 

Since only relative phases are important for oscillations, we can subtract the quantity 
{l/4:E){ml + m|) + {l/y/2)GFNe from the diagonal, and the evolution equation becomes: 



d 



dx K^b 



-A{x) B \fUa 
B A{x) )\vb 



with 



A{x)^^cos{29)-^Ng{x) 



(7.3) 



(7.4) 



For the case of constant density this leads to an oscillation probabihty 

P(i/a Ub) = sin^(26'„) sin^(a;L) 

where 



UJ 



sm'^{29) + cos(2^) 



(7.5) 



(7.6) 



(7.7) 



gives the effective squared mass difference, divided by 4:E, in matter, and 9^ is the relevant 
effective mixing angle in matter, specified by 



sin2(2^„ 



sin2(2^) 



sin^(2e) + (cos(2e)-^^^^^)' 



(7. 



Thus, the resonance condition is 



E^13 GeV 



Am^ 



3 X 10-3 eV^ 



3 g/cm^^ fl/2 



cos29 



(7.9) 



where we have introduced scaling factors normalized by typical values of the density and the 
fraction Yf. = Z/A in the upper mantle. 

Letting the vacuum oscillation length L^ac be defined as Ly^c = 47r£^/|Am^|, the effective 
oscillation length in matter, defined by uiL^ = vr, is 



sin2(2^) + cos(2^) 



2^/2GFNeE\'^ 
Am^ 



'1/2 



(7.10) 
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We recall that, as is evident from these formulas, the oscillation probability in matter depends 
on the sign of cos 20, i.e., whether 9 is in the first or second octant, given that one takes 
Am^ > 0. 

We next recall the formulas for matter effects on oscillation probabilities in the three- 
flavor case with the one-Am^ dominance approximation. Here, the evolution of the weak 
eigenstatcs is given by 



where 



.d_ 
dx 



(7.11) 









"'I 










Us/ 



(7.12) 





V 






2 

mi 



(7.13) 



V2GFNe 
V= \ 




(7.14) 



Subtracting ml from the diagonal, becomes 



/O \ 

Am'i^i 
VO Aml^ + AmlJ 



(7.15) 



In order to calculate the oscillation probabilities for long-baseline terrestrial neutrino oscilla- 
tion experiments and for analysis of atmospheric neutrino data, it is convenient to transform 
to a new basis defined by (e.g. [34]) 



u = R23I' 



(7.16) 



The evolution of is given by 



1 

2E' 



In the one-Am^ approximation, this can be reduced to 



( ^sl^Aml^ + ^/2GFN^ ^si3Ci3Ami2e-^^\ 



H ^ 







2E 



'e " 2E' 










^ siaCisAmlne''' ^c^^Aml 



32 



(7.17) 



(7.18) 
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It can be seen now that in the basis {ug, u^, Pt) the three-flavor evolution equation decouples, 
and it is enough to treat the two-flavor case. We deflne S and P by 




(7.19) 



and 



P = 




(7.20) 



Transforming back to the flavor basis (z/g, ^t)i the probabilities of oscillation become 



If in (7.18) we subtract from the diagonal the quantity D = (l/4E)Ami2 + {l/V2)GFNe, 
we see that it is then necessary to solve the evolution equation for a two-flavor neutrino 
system as in equation (2.11), where in A and B, Am^ = Am32 and 9 = 9i^. For the case of 
constant density, S = e"*^^, so that 5*33 = e~*^^(cosa;L — i{A/u!) sina;L) and P is given by 
equations (7.6)-(7.8). Explicitly for t/^ — > Ug, 



In this case, as in the two-flavor analysis, the interaction with matter makes the oscilla- 
tions sensitive to the sign of Am^^^. For antineutrinos, the matter potential has the same 
magnitude and opposite sign, so one has to solve the same evolution equation where V is 
replaced by —V. Consequently, if the oscillation probabilities are enhanced by the presence 
of the matter for neutrinos, as they are for Am^^^ > 0, then they will be suppressed for 
antineutrinos and vice versa. From these results it is evident that V ^ —V is equivalent 
to Am^ — > — Am^. The neutrino factory physics program intends to use this property to 
obtain the sign of Am32 by comparing P(fe — > 1^^) and P(Pe — ^ i^n) (e.g., [30]). 

Let us now present some numerical results for matter effects in the case of three-flavor 
oscillations with two Am^ values. We concentrate on the case of large pathlengths and the 
J^ij. transition relevant to the existing data on atmospheric oscillations. For simplicity, 
we take the CP- violating phase equal to zero here, but it is straightforward to include it (see 
below). These results can also be applied to data on Ug — > that might become available 
with a possible future neutrino factory. In Fig. 1 we plot Pii^fi I'e) as a function oi L/E 
for sin^ 2^23 = 1, Am|2 = 3 x IQ-^ eV^, sin^ 2^13 = 0.04, sin^ 2^12 = 0.8, and Am^^ = 2 x 10"^ 
eV^, the upper end of the LMA region. The higher- frequency oscillations are driven by the 
terms involving sin^ ^32 while the lower-frequency oscillation is driven by the terms involving 
Am|^. The one-Am^ approximation is shown as the dashed curve; of course, this lacks the 
low-frequency oscillation component. One sees that the full calculation differs strikingly from 
the result of the one-Am^ approximation. 




(7.21) 
(7.22) 
(7.23) 



P(i/^ i/e) = Pii^e T^n) = sin^(26li3,„i) sin^(6l23) sin^(a;32,m-^) 



(7.24) 
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FIG. 1. Plot of P(i/^ 



as a function of L/E for sin^ 6*23 = 1, Am§2 = 3 x 10"^ eV^ sin^ 26*12 = 0.8, 



and Aml^ = 2 x IQ-^ eV^, and sin^ 29i3 = 0.04. 

Even for the best-fit LMA solution, the effect of Am^i can be large for large pathlengths, 
and this would affect the z/g oscillations in atmospheric neutrino data, as shown in 

Fig. 2, for which we take the central values of sin 26^21 and Aml^ in the LMA fit, (3.8) and 
other parameters the same as in the previous figure. Note that for the dominant z/^ Vr 
transition in the atmospheric neutrinos, Am^i effects are not so important; this is clear from 
the fact that this transition does not directly involve z/g. 

0.5 



full result, Am 21=5x10 eV 
one mass scale 



0.4 



^^0.3 

> 
t 

Q- 0.2 



sin (2e|>0.8 
sin'(2e,3)=0.04 



0.1 



0.0 




5000 10000 

L /E[km/GeV] 



15000 



FIG. 2. Plot of P{v^ -f Ve) as a function of L/E for sin^ 6*23 = 1, Am§2 = 3 x 10"^ gya^ central LMA 
values sin^ 26*12 = 0.8 and Am^^ = 5 x 10^^ eV^, and sin^ 26*13 = 0.04. 
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We next show, in Fig. 3, the result of integrating (7.11) in the full three-flavor mixing 
scenario and using the actual density profile of the Earth as given in [35]. For this figure 
we use sm^2e2z = 1, ^ml^ = 3 x lO'^ eV^, sin2(2^i2) = 0.8, Am^^ = 5 x 10"^ eV^, and 
sin^(26'i3) = 0.04. As expected, the Am^i corrections are big for low energies and large 
distances. For the choice of a large distance, L = 10*^ km (shown in fig. 3), we observe a 
very significant difference between the full calculation and the one-Am^ approximation. This 
shows again (as does the recent illustrative study in Ref. [18]), that it would be valuable to 
carry out a more complete analysis of the SuperK and other atmospheric neutrino data with 
not just three-flavor oscillations, but also two Am^ values included. Although the SuperK 
fit to its data shows that the z/g oscillations make a small contribution, it is important 

to include this contribution correctly, and the one-Am^ approximation is not, in general, 
reliable for this transition. 



full result, Am\ =5xlO"'*eV^ 
one mass scale 



L=10000Km 

sln'(2e,3)=0.04 




0.00 ^ 



E [GeV] 



FIG. 3. Plot of P{i^f, i/e) as a function of E for sin^ 6'23 = 1, ^mj^ = 3 x lO'^ eV^, and L = lO'^Km, 
with other input values as shown. This calculation takes account of the full density profile of the earth. 



VIII. CP VIOLATION 

It is well known that there can be substantial leptonic CP violation observable in neutrino 
oscillations, in the context of three-flavor mixing and that this depends on both Amgg and 
Am^i being nonzero. The observation of leptonic CP violation is, indeed, a major goal of a 
neutrino factory [30]- [38]. Our main point here pertains to the accuracy with which input 
parameters might be known at a time when a neutrino factory might operate, which, in turn, 
leads one back to the necessity of a general three-flavor, two-Am^ analysis of atmospheric 
neutrino data and data from studies of u^j, —>■ oscillations with an intense conventional 
neutrino beam. To elaborate on this, we recall that at a neutrino factory, a potentially 
promising way to measure CP violation is via the asymmetry 
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This method has two main appeals: (i) one can produce equally intense initial fluxes of 
and by switching the stored beam between and \x^\ and (ii) a detector for long-baseline 
neutrino oscillation searches with a neutrino factory will have the ability to identify outgoing 
//^'s and measure their electric charges, and its detection efficiency will be equal for the two 
signs, so no bias will be introduced in this measurement. The complication with this method 
is that, even in the absence of any intrinsic CP violation, the asymmetry does not vanish 
because matter effects reverse sign between neutrino and antineutrinos, and the earth is 
not CP-symmetric. Therefore, the challenge with this method will be to determine these 
matter effects with sufficient accuracy to be able to disentangle them from the intrinsic CP 
violation. An important source of information here will be the anticipated measurement 
of ^^13 from the JHF-SuperK experiment on z/g oscillations, since matter effects are 

sensitively dependent on this parameter [33]. As we have shown in a previous section, for 
an accurate determination of ^13, the one-Am^ approximation is not, in general, reliable, 
especially if sin^ 2^21 and Am|i are near to their maximal values in the LMA solution to the 
solar neutrino deficit. 

It should be noted that there are also plans to try to measure CP violation with an 
intense conventional beam, e.g. in the JHF-SuperK experimental program, by comparing 
overall rates of Vf. and — > Pe and also energy dependences of these signals [26,27] (see 
also [38]). For sufficiently large sin^ 2^13, sin^ 2^12, and Am|i, these methods might provide a 
way to measure CP violation complementary to that used with a neutrino factory. There are 
a number of challenges with this approach: (i) the fluxes of z/^ and are different; (ii) the 
event rates would involve several different cross sections, Vf.n — > e~p and PgP e"'"n in oxygen 
nuclei, as well as v^p — > e^n on the hydrogen nuclei in the water molecules; (iii) since it is 
not possible to determine the sign of the with SuperK, the comparison would have to be 
done with data from different periods of operation. Assuming these experimental challenges 
can be met, the importance of accurate inputs for the various parameters is evident from 
the formula for the CP- violating asymmetry 

A, = ^S'^^ ^ - ^ ^ (8.2) 

For the JHF-SuperK basehne matter effects are small, so we shall consider the expression 
for this asymmetry in vacuum (where A2 = —Ai). Using 

-P(i^M ^e) - P{i^ix i^e) = 4J(sin2032 + sin 2021 + sin 2^13) 

= 16 J sin 032 sin ^31 sin 02i (8.3) 

where J was given in eq. (2.3), substituting the expression for P{i'n fe) + P{^ix ~^ ^e) 
from (4.1), and using Am|i << Am|2, one has 

sin(26li2) cot 6'23 sin 5 sin 02i . . 

A2 ^ — ^ — ^-7 8.4 

sm 6113 
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Thus, to extract an accurate measurement of sin S, it is clearly important to have sufficiently 
accurate inputs for quantities such as ^13, and this, in turn, motivates the generalization 
of the one-Am^ approximation that we have presented above for the analysis of data on 
neutrino oscillations. 

IX. CONCLUSIONS 

In this paper we have performed calculations of neutrino oscillation probabilities in a 
three-flavor context, taking into account both Am^^^ and Amg^i scales. We have shown that 
for values of sin^(2^i3) ~ 10~^ in the range of interest for long-baseline neutrino oscillation 
experiments with intense conventional neutrino beams such as JHF-SuperK and with a 
possible future neutrino factory, and for Am^^j ~ 10""^ cV^, the contributions to z/^ 
oscillations from both CP-conserving and CP- violating terms involving sm"^ {AmlgiL / (AE)) 
can be comparable to the terms involving sin^(Am^^^L/(4ii^)) retained in the one-Am^ 
approximation. Accordingly, we have emphasized the importance of performing a full three- 
flavor, two-Am^ analysis of the data on — > z/g oscillations from an experiment with a 
conventional beam, and on z/g — z/^, z/g — > z/^ oscillations from experiments with a neutrino 
factory. In our study we have included calculations of matter effects in a three-flavor, two 
Am^ framework. Our results also motivate the analysis of atmospheric neutrino data in this 
generalized framework. 
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X. APPENDIX: ANALYTIC APPROXIMATION FOR THREE-FLAVOR 
TWO-Am^ OSCILLATIONS IN MATTER 

In matter, for constant density, it is still possible to solve the evolution equation exactly 
( [36]). However, these results are rather complicated. For fits to the data or other studies, 
it can be useful in this case to use approximate formulas that can still very well describe the 
oscillations. 

For the LMA solution to the solar data, the effects of Aml^i can no longer be neglected, 
but there are cases where they are small in long-baseline experiments for which the dominant 
oscillation is controlled by Am^^^^. In these cases, one can thus treat the effects of Am^^^ as 
a small perturbation. This has been done in [37], where the matter effect is also take to be 
a small perturbation. This is a good approximation at short and medium distances. Here 
we give generalized formulas that treat the Am^^; as a perturbation but allow large matter 
effects, as is necessary in very long baseline experiments. 

In order to calculate oscillation probabilities it is now convenient to work in a basis 
defined by: 
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The evolution of v' is given by 
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From this we obtain ^'{x) — 5'V(0), with 
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The first term gives the one-Am^ contribution, which, after the rotation to the flavor 
basis (i/g, f^, Vr), is the same as the one obtained in the previous section. The second term 
contains the corrections due to the Am^^i term. 5"^ is given by 
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After we rotate back in the actual flavor basis, the oscillation probabilities will be given by 

Pii^a - lyb) = mR\,S'RuU^)ab\' (10.15) 

The results from Eq. (10.15) are compared in Fig. 4, 5 with those obtained from the 
exact numerical result and with the result for the one-Am^ approximation for the following 
set of input values: sin^ 2^^23 = 1, Am^g = 3 x 10~^ eV^, sin^ 2^^soi = sin^ 26'i2 = .8 and 
^^"soi — '^^21 — 10~^eV^ (or Am^^i = for comparison), and sin^ 26*13 = 0.04. 
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FIG. 4. Plot of P{i'e v^) as a function of E for L = 2900 km. Other input parameter values are 
given in the text. 
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E [GeV] 

FIG. 5. Plot of P(z^e v^) as a function of E for L = 2900 km, with detail of the region from i? = to 
E = 1 GeV. Input parameter values are given in the text. 
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